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A high-energy spin resonance mode is known to exist in many high-temperature superconduc- 
tors. Motivated by recent scanning tunneling microscopy (STM) experiments in superconducting 
Bi2Sr2CaCu20g+5, we study the effects of this resonance mode on the local density of states (LDOS). 
The coupling between the electrons in a d-wave superconductor and the resonance mode produces 
high-energy peaks in the LDOS, which displays a two-unit-cell periodic modulation around a non- 
magnetic impurity. This suggests a new means to not only detect the dynamical spin collective 
mode but also study its coupling to electronic excitations. 
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A prominent feature in the excitation spectrum of the 
high- Tc superconductors is the "41 meV" collective spin 
resonance mode, seen by inelastic neutron scatteriiig ex- 
[ periments in most of the cuprate families P, 0, Bi 0- 
The physics of this resonance mode - including its mi- 
croscopic origin, its connections with other physical prop- 
erties, as well as its role on superconductivity itself - 
has been the subject of considerable debate. Given the 
recent developments of the atomic resolution scanning 

tunneling microscopy (STM) SHU 1113, it is timely 

to address the possible manifestation of this resonance 
mode in the local densit y o f states (LDOS). A number 
of theoretical works 0, llll lisj have addressed the 
effect of related spin fluctuations on the LDOS. These 
works focused on the pinning of the spin fluctuations by 
impurity: a dynamical spin mode centered around the 
wavevector Q leads to a 2Q spatial modulation in the 
low-energy LDOS. This result is smoothly connected to 
what happens in the case of a static spin-density-wave or- 
dering |l4ill5l| . However, for the resonance mode - which 
is sharply peaked at Q = (tt, tt) - such effects would not 
be manifested [since 2Q = (27r, 27r) is equivalent to (0, 0)]. 
There are also works about quasiparticle scattering from 
nomagnetic impurities |l6lll7l |. 

In this paper, we show that the coupling of d-wave 
quasiparticles to the resonance mode does produce spa- 
tial modulations in the LDOS around an impurity. The 
feature is located at relatively high energies, « ±(Ao -t- 
f2o)j where Aq is the maximum superconducting en- 
ergy gap and CIq the resonance energy. In addition, 
the wavevector of the LDOS modulation is close to 
Q = (tt, tt). Our predictions could be observable by 
the STM experiments. Such STM studies represent a 
new means to characterize the coupling between the elec- 
tronic excitations and the resonance mode. The STM 
feature we discuss relates to the "peak-dip-hump" struc- 
ture of the angle resolved photoemission spectroscopy 
(ARPES) m El m lii, il; the inference about 
the electron-spin coupling from the ARPES and related 
spectroscopies is a topic of recent controversy 11101 and 



we hope that the STM studies we propose will shed new 
light on this important issue. 

We start with a model Hamiltonian describing two- 
dimensional electrons coupled to a collective spin mode 
and in the presence of a single-site impurity: Ti. = 



"Hbcs + 
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n 



imp • 



Here the BCS-type Haixiiltonian 
for a uniform d-wave superconductor is given by TLbcs = 

where c\^^ (cko-) creates (annihilates) a conduction elec- 
tron of spin a and wavevector k, is the normal state 
energy dispersion for the conduction electrons, ^ the 
chemical potential, and Ak = -^(cosA:^ — cosfcy) the d- 
wavc superconducting energy gap. The coupling between 
the electrons and the resonance mode is modeled by an 
interaction term Hgp ~ ■ ^i, where the quanti- 

ties g, Si, and Si are the coupling strength, the electron 
spin operator at site i, and the operator for the collec- 
tive spin degrees of freedom, respectively. The dynam- 
ics of the collective mode will be specified below. The 
impurity scattering is given by H^mp = UoY^^cl^coa, 
where without loss of generality we have taken a single- 
site impurity of strength Uq located at the origin, = 0. 
By introducing a two-component Nambu spinor oper- 
ator, ~ (ckti c1ij.|)"^; the matrix Green's function 
for the d-wave BCS Hamiltonian "Hbcs is determined 
i^n - $k -Ak 

-Ak iujn + Ck 

Sk — fi and ujn = i'2n + 1)ttT is the fermionic Matsubra 
frequency. We have also assumed that the d-wave pair 
potential is real. 

For a homogeneous system, where only the inelastic 
scattering of quasiparticles from the collective mode oc- 
curs, we calculate the self-energy to the second order in 
the coupling constant (see Fig. ^): 



by Gg \k;iw„) 



, where ^k 



E(k; iuJn) 



X(q; i^i)Goi'k - q; iuj,, - i^i) , 

q f^; 



where x(q;*^i) is the dynamical spin susceptibility 

Xy(T) = (TV(5'f (r)5J(0))) and Qi = 2/7rr the bosonic 
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FIG. 1: a) Diagrams for G(p,p + q;aj). The thick sohd hne 
represents the full conduction electron Green's function in the 
absence of impurity, Gg. The dashed line is the coupling to 
an impurity, specified by a cross. The wavy line denotes the 
propagator of the collective mode; b) Self-energy diagrams for 
the conduction electrons in the absence of impurity, E. The 
thin solid line is the bare (BCS) conduction electron Green's 
function. Go. 



Matsubra frequency. The dressed Green's function is: 



(1) 

The corresponding real-space dressed Green's function 
Gj){i,j; iuJn) is obtained through a Fourier transform with 
respect to — Vj. For the d-wave pairing symmetry, 
one can show that the local Green's function, g^{iu!n) = 
GQ{i,i\iu!n) is diagonal. In the presence of a single-site 
impurity at = with potential strength [/qj the site- 
dependent Green's function can be written in terms of 
the T-matrix: 

G{i.3;E) = a,{i,r,E) 

+ Y.Go{i,l;E)fimiE)Go{m,j;E) .{2) 

Im 

Due to the vertex corrections induced by the coupling to 
the collective modes (Fig. QJi), the T matrix in general 
contains site-off-diagonal terms. We will first carry out 
the calculation without the vertex corrections, in which 
case fim = fdi^o6jn,o, with f-^ = UQ^a^ - g^, where 0-3 
is the 2-component of the Pauli matrix. The LDOS at 
the i-th site, summed over two spin components, is 



Pin,E) 



-ImGii(i, i;E + i^) , 



(3) 



where 7 = 0^. 

Up to now, our discussion and formulation are quite 
general and can be used to study the effects of any dy- 
namic mode once the susceptibility x is known. We treat 
the susceptibility in a phenomenological form (based 
on the inelastic neutron scattering observations), see 
also El: 



x(q;i^z) 



1 



1 



iVli — Via i^i + f^o 



(4) 



where we denote the wavevector Q = (7r,7r) and the 
mode energy by fig. This form is especially suitable for 
the optimally doped YBa2Cu306+y (YBCO) compounds 
in the superconducting phase, where the observed neu- 
tron resonance peak is almost resolution-limited in en- 
ergy and fairly sharp in wavevector. The resonance peak 
in BSCCO is broadened in both energy and wavevector. 
The finite width in the wavevector space might be impor- 
tant for the ARPES lineshape in general and in particular 
the understanding of the ARPES spectra away from the 
M points [k = (tt, 0) and symmetry-related points] of 
the Brillouin zone |22| , but should not change the quali- 
tative conclusion of our work: the LDOS effects we will 
discuss arise from the fact that the dominant effects of the 
resonance mode on the single-electron spectral functions 
occur near the M points which is expected to remain to 
be the case beyond our simplified form for the suscep- 
tibility. In addition, given that the peak in BSCCO is 
still quite sharp in energy, we expect that the main effect 
of the broadening in energy of the resonance mode is to 
extend the bias window for the LDOS feature we will dis- 
cuss. We have also neglected the incommensurate peaks 
seen in the inelastic neutron scattering experiments in 
YBCO (the part that disperses "downward" away from 
the resonance peak) |2J, VZA I2(t[ l27l | , since their spectral 
weight is significantly smaller than that of the resonance 
mode. For the normal-state energy dispersion, we use 
ek = — 2t(cos kx + cos ky) — At' cos kx cos ky, where t and 
t' are the nearest and next-nearest neighbor hopping inte- 
gral. Unless specified explicitly, the energy is measured in 
units of t. We choose t' = —0.2 to model the band struc- 
ture of the hole-doped cuprates. Since the maximum 
energy gap for most of the cuprate superconductors at 
the optimal doping is about 30 meV while the resonance 
mode energy is in the range between 35 and 47 meV, we 
take Aq 0.1 and Qq = 0.15 (i.e., 1.5Ao). The chem- 
ical potential is tuned to give an optimal doping value 
0.16. To mimic the intrinsic life time broadening, in our 
numerical calculation we take 7 of Eq. Q to be 0.04Ao. 
A system size of x Ny ^ 1000 x 1000 is taken in the 
numerical calculation. 

In the absence of impurities, the density of states 
is the summation of the spectral function, Ak{E) = 
— -llmGg ii(k; E + 17), over all wavevectors k. Fig. El^a) 
shows the spectral function at an M point of the Bril- 
louin zone. Without the electron-mode coupling, as 
is well known, the spectral function is peaked at the 
maximum gap edges ±Ao. As the electron-mode cou- 
pling is switched on, new peaks emerge at the energies 
±Ei w ±(Ao + r^o)- (For simphcity, we have neglected 
the broad "background" part of the single-electron spec- 
tral function in our consideration.) The peaks in orig- 
inate from the poles of the Green's function Gq h- Note 
that the weight of the peak at — Aq is larger than at 
Aq because the van Hove singularity is below the Fermi 
energy. Since the spectral weight of the spin resonance 
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FIG. 2: Spectral function Ai^{E) at the wavevector k = (tt, 0) 
for g/Ao — 3/\/2 (a). The density of states, X]k^k(-E'), is 
shown in (b), for various values of g/Ao = 0, 1/V2, 2/^2, 
and 3/V2 (from lower to upper); here, for easier viewing, the 
consecutive curves are shifted by 0.2 along the vertical axis. 



mode [i.e.; Imx(q;tj)] is peaked at Q, the feature of the 
quasiparticle self-energy is the strongest around the M 
points of the zone because they are connected by Q. The 
singularity in the quasiparticle self-energy causes addi- 
tional poles in the Green's function. As the coupling 
constant g increases, these peaks are shifted to higher en- 
ergies and, in addition, their spectral weight is enhanced; 
simultaneously, the weight of the superconducting coher- 
ent peaks is reduced to obey the sum rule. The shift 
of states due to inelastic scattering is expected in DOS 
and is also expected for scattering off local mode |28| . 
Fig.l^b) plots the density of states and clearly shows that 
the high energy peaks still occur around ±Ei. In other 
words, the contributions from near the M points domi- 
nate the wavevector summation for the density of states, 
reflecting the flat nature of the normal state band near 
this point. Furthermore, the highly asymmetrical struc- 
ture in the DOS at energies —Ei and Ei comes from the 
singular structure in the quasiparticle self-energy. These 
results, for the clean case, are consistent with earlier stud- 
ies of the ARPES [H El 113 El El and DOS E^. 

We are now in a position to address the LDOS in the 
presence of a single nonmagnetic impurity. For concrete- 
ness, we take the on-site potential Uo — IOOAq. Fig. |2| 
shows the LDOS directly at the impurity site, as well as 
at its nearest neighbor. The near-zero energy resonant 
state trig gere d by the quasiparticle scattering from the 
impurity |22j is robust against the electron-mode cou- 
pling. Our key new results are two-fold. First, the impu- 
rity modifies the shape of the spectral features at ±Ei, 
which can now be either a dip or a peak. Second, and 
more importantly, these high energy features at ±Ei ex- 
hibit a spatial dependence. At the impurity site, the 
LDOS displays a dip at —Ei but a peak at +Ei. The 
behavior is reversed at the site closest to the impurity. 

To explore this spatial variation of the LDOS in more 
detail, we have calculated the LDOS, with the energy 
fixed at —El, in the vicinity of the impurity with and 
without the mode coupling. In the absence of mode cou- 



FIG. 3: Local density of states at the impurity site (a) and 
at its nearest neighbor (b). The coupling constant g/Ao = 
3/^/2. 




FIG. 4: The spatial variation of the LDOS around the impu- 
rity with the mode coupling g/Ao = S/y/2 at E = —E\ (a), 
and its corresponding Fourier spectrum p(q, — iSi) (b). The 
density of states is measured with respect to its spatial aver- 
age value. The image window of 25 x 25 plaquettes is taken 
from a system of size 1000 x 1000. Here the impurity-induced 
Friedel oscillation along the diagonals have been filtered away; 
see main text for details. 

pling, we obtain the results (not shown) which are con- 
sistent with previous studies E3- the LDOS exhibits a 
Friedel oscillation along the diagonals of Cu02 plane but 
no other non-trivial features. When the mode coupling 
is turned on, in addition to the Friedel oscillation along 
the diagonals, the LDOS displays a new type of modu- 
lation with a period of 2a in the wide regions along the 
bond direction. This new modulation can be more easily 
seen Fig. Efa) when the pre-dominant Friedel oscillation 
is filtered away. It occurs in four disconnected triangles 
in the field of view. 

Also in order to highlight the new modulation, we 
find it useful to perform a filtered Fourier transforixi, 
p{t\,E) = Y^^e^^'^^ p{vi,E), where denotes a sum- 
mation over all the sites in four triangles in Fig. |^a). 
The resulting Fourier-transformed image E3| is given in 
Fig. EJb) , which unambiguously shows that the new fea- 
ture induced by the coupling to the spin resonance mode 
has a spatial modulation wavevector (tt, tt). 

This new type of LDOS modulation with a wavevector 
close to Q = (tt, tt) reflects the dominance of the collec- 
tive mode effect near the M points of the Brillouin zone. 
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It follows from Eqs. H2I3() . with the form of the T-matrix 
under consideration and when the chosen field of view 
has an inversion symmetry with respect to the impurity 
site, that the Fourier transformed LDOS for any finite q 



IS, 



-Im / lip 



Go(p + q;i?)f(ii;)Go(p,i?) .(5) 
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At positive energies, the poles shown in Fig. |2] for 
Gq{-9\E) are located [H El El at and 

-Bp-Q + ilo, respectively. (Here, = \JQ, + A^.) Like- 
wise, Gq(p + q; E) is the sum of two poles, one at -Ep+q 
and the other at i?p+q_Q + Oq. At = i^i, Eq. © 
is dominated by the convolution between the poles at 
Ep-Q + rio and Ep+q-Q + flo. This term in turn is dom- 
inated by the contributions corresponding to when both 
p — Q and p -I- q — Q are near to the M points, leading 
to a /5(q, —El) that is peaked near {tt,7t). 

Similar phase-phase considerations show that the ver- 
tex correction terms lead to a similar momentum depen- 
dence. In the presence of vertex corrections, the T-matrix 
satisfies an integral equation in the wavevector space. 
The vertex correction to the T-matrix to the same or- 
der (g^) of our calculation for the self-energy is shown in 
Fig. ^a). It involves a wavevector convolution of a form 
similar to that given in Eq. We therefore expect [sJl 
an additive contribution to /9(q, —Ei) that is also peaked 
near (tt, tt). 

Finally, we remark on issues which go beyond the ideal- 
ized model we have considered so far: (i) In the presence 
of a gap inhomogeneity at the nanoscale, the mode signa- 
ture would appear at (Aq) -I- wq, with a slight smearing. 
Here, (Aq) is the spatially averaged superconducting gap. 
The fact that a well-defined peak-dip-hump structure ap- 
pears in the break-junction tunneling spectrum im- 
plies that the smearing is not too large; (ii) The mode 
signature depends on the detailed band structure. In the 
cupratcs, the band is flat near the antinodal points (i.e., 
(tt, 0) etc.), and a mode with a momentum significantly 
different from (tt, tt) will have a weaker effect compared 
to that of the (tt, tt) mode we addressed; (iii) Tunneling 
matrix elements need to be taken into account in order to 
understand the detailed spatial variation of the LDOS as 
observed around zinc impurities in BSCCO . Such 

a filtering effect, however, will not affect our conclusion 
on the momentum of the LDOS modulation. 

To summarize, we have studied the effects of the mag- 
netic resonance mode on the tunneling spectrum in the 
presence of a nonmagnetic impurity. The LDOS around 
the impurity displays resonant features at relatively high 
energies [close to ±(r2o + Ao)], which modulates in space 
with a wavevector close to (7r,7r). Our prediction can be 
tested straightforwardly by operating the existing high 
resolution STM at a relatively high energy window. Such 
studies should shed considerable new light on the physics 



of the spin resonance mode, in particular its coupling to 
the electronic excitations. 
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